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Abstract 

We investigate theoretically the thermal escape behavior of trapped magnetic fluxons in long an- 
nular Josephson junctions in dc magnetic fields, and perturbed by a probing ac current. The study 
is motivated by recently published experimental data that show multi-peaked escape distributions 
for increasing bias current in the extreme low temperature regime when the system is perturbed by 
an ac current. We demonstrate that the observed behavior of multi-peaked escape distributions can 
be reproduced and predicted in the entirely classical, thermally driven sine-Gordon model, which 
is widely accepted as accurately describing the experimental system. We interpret the observed 
multi-peaked distributions as being directly induced by dynamical resonances between the applied 
ac perturbation and the natural oscillation frequency of a trapped fluxon. 

PACS numbers: 85.25.-j, 85.25.Cp, 74.50.+r 
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Several experiments have suggested^ that Josephson junctions, operated in their zero- 
voltage state and close to the critical current, will reveal signatures of quantum levels if 
perturbed by ac currents at frequencies which, when multiplied with Planck's constant, cor- 
respond to the energy between the anticipated levels. The signatures are found by slowly 
increasing the dc bias current of the junction and recording the current at which the voltage 
across the junction switches from zero to a non-zero value. By repeating such experiments, 
one can obtain a distribution of critical switching currents, which reflects aspects of the 
nature of the switching device. These studies have been inspired by extensively conducted 
work over the past three decades confirming the consistency between Kramers theory for 
thermal escape with the predicted behavior of superconducting devices at non-zero tem- 
perature for various kinds of Josephson junction systems ^ 4 ! 5 ' 6 ! 7 ' 8 . Given the consistency in 
the classical, thermal regime of Josephson junctions, similar escape measurements in the 
low temperature regime, where quantum uncertainty should dominate thermal fluctuations, 
have been conducted. These have provided the possibility for non-thermal features in the 
measurable switching distribution. A specific aim has been to identify signatures of quan- 
tum levels expected in the extreme low temperature regime of the zero-voltage junction. 
Such energy levels could manifest themselves through the escape statistics if more than one 
energy level is statistically possible, since a junction in an excited energy state will have a 
higher probability of switching at lower bias currents than a junction in the ground state. 
Thus, if multiple quantum levels can be populated, then the switching current distribution 
may reflect this as a multi-peaked profile, each peak identifying a unique quantum level. 
Experimental efforts have identified and reported such multi-peaked switching distributions 
in a variety of Josephson junction systems with application of a small ac current. The ra- 
tionale for this ac perturbation is that it provides the energy (at a particular frequency) to 
populate specific excited energy levels. 

A particularly noticeable set of observations was reported for an annular long Josephson 
junction in a static magnetic field, and with a single trapped magnetic flux quantum^. This 
system is know n 1 ^ 11 to exhibit behavior similar to that of a small-area Josephson junction, 
and thus, should behave equivalently to small area junctions when conducting switching 
measurements 6 . Since a magnetic fluxon is a macroscopic object, anomalous switching 
behavior related to the fluxon dynamics would imply that macroscopic quantum dynamics 
had been identified. It is the aim of this paper to investigate what can be expected from an 
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equivalent purely classical model of the system, and compare the results to the published 
experimental results to assess the interpretation of the experiments. Thus, the physical 
object of our investigations is an annular Josephson junction of inner radius R and width W 
in the limit that the length of the junction 2ir(R + W/2) is much larger than the Josephson 
penetration depth^ . The total area of the junction is (irW 2 + 2ttRW) and therefore, if the 
width is much less than the inner radius, a good approximation for the area of the annular 
junction is 2nRW. We make this last assumption, but we also assume that the width of 
the junction is smaller than the Josephson penetration depth, meaning that we restrict our 
interest to one- dimensional physics along the barrier forming the junction. In this limit, 
relevant electrical parameters are the capacitance (C), inductance (£), and the sub-gap 
conductance (Q) per unit length of the junction. For the case of the annular junction that 
we consider, C = eos r W/T and C = ^qD/W. Here, e$ and /xo are, respectively, the vacuum 
permittivity and permeability and e r is the relative dielectric constant of the oxide of the 
junction. The parameter T represents the thickness of the oxide barrier while the parameter 
D = T + Ai + A2 is the magnetic thickness of the junction. These parameters are all linked 
in Josephson junction through the electromagnetic wave velocity l/y/CC = a/ (T /V) / n eoE r 
in the oxide barrier. 

We study the perturbed sine-Gordon model in the following form, 



where <p(x, t) is the phase difference between the quantum mechanical wave functions of the 
two superconductors defining the junction. Equations (H}-(HJ) model an annular Josephson 
junction of inner radius R and width W when a dc magnetic field is applied in the plane 
of the junction and the forcing terms include dc and ac currents as well as thermal noise. 
Assuming uniform Josephson supercurrent pair density j c over the area of the junction, 
the critical current density (per unit length) is Wj c . This is also the characteristic current 
density against which all current densities are measured in the normalized equations (HJ-fJU). 
Tunneling of Cooper pairs of electrons is described by the simp term length. Dissipative 
quasi-particle current density is represented by the term proportional to a = Qhu p /2eWj c 



<Ptt - <Pxx + sin Lp 





<p(x,t) 



Vdc + Vac sin ujdb + Tk sin kx — a<p t + n(x, t) 



ip(x + l,t) + 2tt , 



(1) 
(2) 
(3) 
(4) 
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where u p is the unperturbed plasma frequency u p = y 2 c and where $o = ^ ~ 2.07 • 
10~ 1B Wb is the flux quantum. The spatial coordinate x is normalized to the Josephson 



penetration depth Aj = y j^cwf w hile the temporal coordinate t has been normalized to 
ujp 1 . An external magnetic field in the plane of the junction (see Fig. is represented 
by r ^ , which is normalized to the characteristic field Hq = j c Aj. Given the annular 
geometry, the effect of the magnetic field is a spatially harmonic potential function with 
periodicity I = 27T ( R +w/2) ^ e normalized length of the junction; thus, k = 2n/l. The 
junction is driven by a dc current, an ac current, and a thermal noise current from the 
dissipative quasi-particle term proportional to a, represented by the normalized quantities 
Vdc, Vac and Ud, as well as n, where n is given by the dissipation fluctuation relationship of 
Eqs. (J2J) and dH^. The normalized temperature 9 is measured relative to the characteristic 
temperature T = Ej/ks, where Ej = $0 2 ^ Jlc is the Josephson energy density. Finally, 
since we consider only the situation of a single trapped flux quantum, we always impose the 
boundary condition of Eq. Q . 

As outlined in Refs . 10 ' 11 , the dynamics of an annular Josephson junction with a single 
trapped fluxon (flux quantum) should behave similarly to a small-area Josephson junction, 
which in turn is typically modeled by a pendulum equation 12 . This can be illustrated by 
defining the wave momentum P, 

P = ~l I VxVtdx , (5) 



and using the standard kink-soliton solution, 

X — Xq — ut 



ip s = 4 tan 1 



exp 



v. 



(6) 



where u = Xq and the resulting fluxon momentum is P = u/y/1 — u 2 . We then obtain the 
following dynamical equation for the fluxon position, xq: 



dP n n . 71 ( 1 i A 

— — = —f]dc H — Vac sin u d t — aP 1 Ksech -Tikyl — u l sin kx 

at A A 4 \2 J 

The thermal fluctuation term associated with the fluxon 



o + v(t) • (7) 



M*)> = o (8) 
(HtHt')} = l-^L=5(t - t') . (9) 



4 



In the following, we restrict the analysis to the non-relativistic limit, y/1 — u 2 — > 1. In this 
limit, the above equation is equivalent to the usual pendulum model for a small Josephson 

junction, and for 6 = a = r] ac = 0, we can directly determine the critical current r) c and 

( f ) 

linear plasma oscillation of the ffuxon ujp . 



where r]c is the maximum dc bias current for which zero-voltage states are possible in the 

( f) 

system, and ujp is the small amplitude oscillation frequency of the fluxon in the washboard 
( f) 

potential for 77^ < ry c . The corresponding equilibrium location of the static fluxon is given 



These expressions are derived for a — n — r} ac — 0, I 1, and small to moderate values of 
T, and arise from a collective coordinate treatment of the fluxon. Thus, the results account 
neither for spatial fluxon deformation due to the perturbation terms nor for the periodic 
boundary conditions. 

The very close analogy between the behavior of a long annular Josephson junction with 
a single trapped fluxon in a magnetic field, and a small area Josephson junction has led 
to several theoretical and experimental studies. One of the latest is the investigation to 
identify macroscopic quantum tunneling, triggered by the application of microwave radia- 
tion, of a flux quantum through the effective washboard potential barrier—. Since it was 
recently demonstrated^ 5 that the signatures in anomalous switching distributions obtained 
from small-area Josephson junctions, generated by the application of ac perturbations and 
attributed to quantum levels, can be reproduced in the classical pendulum model and in 
experiments conducted at temperatures well above the quantum transition temperature, it 
is logical to anticipate the same overlap in signatures between the published experimental 
low temperature results and the classical fluxon behavior. 

To investigate the relationship between the classical model and the experimental observa- 
tions reported in Ref.^, we have conducted extensive numerical simulations of the statistics 




(10) 
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of thermal fluxon escape from the magnetic-field-induced potential barrier in a manner simi- 
lar to the experiments. The system is initiated with a trapped fluxon, a given magnetic field 
T, and a dc bias current well below the critical value, rjc . Thermal noise and an ac bias 
current at a given frequency are also applied. The dc bias current is continuously increased, 
typically with a normalized rate of r)d c < 10~ 6 . The system is defined to have switched from 
the "static" zero- voltage state by meeting three conditions: 1) the fluxon position, measured 
as the position of the largest value of \<p x \, exceeds 1/2, 2) the phase {<p) x — sin -1 rjd c exceeds 
Ti, and 3) the normalized voltage (tp) x exceeds | of the normalized voltage corresponding to 
the perturbation estimate of the steady fluxon motion when V = 0^: 

= ; = • (13) 

r= Vac =n=o v /(4«)2 + (rr Vdc y 

All three switching criteria yield nearly identical switching distributions for r)d c 3.5 • 10~ 5 . 
For the very fast sweeps, one of the criteria may be activated immediately due to transient 
effects from the initial conditions. In such a case, we define switching to have happened 
when all three criteria are fulfilled. The simulations are carried out with a spatial central 
difference discretization of dx = 0.025 and a Verlet temporal integrator with time step, 
dt ~ 0.02. The normalized damping parameter has been chosen here as a = 0.001 for all 
data shown. 

We first show the results of a set of simulations of 75,000 switching events for I = 10, 
T = 1, and Ud = 0.2. Figure 121 shows the resulting probability density p(rjd c ) for switching 
at a given value of rjdc f° r different values of the ac amplitude rj ac . The density is normalized 

such that J pdrjdc = 1. For very small values of rj ac , we find a single-peaked localized 

( f) 

switching distribution for rjdc values just below the predicted critical switching current rj c . 
However, for increasing ac amplitude, we observe a nontrivial transition to more complicated 
statistics of the switching events. For rj ac in the interval of 5 • 10~ 4 < rj ac < 10~ 3 , we 
observe clear evidence of a double-peaked distribution, closely resembling the experimental 
observations of Ref.«, and for ac amplitudes higher than 10~ 3 , we again find a localized 
single-peak distribution, but centered at a lower rjd c value than for rj ac = 0. The origin of the 
multi-peaked distributions observed experimentally has been interpreted as evidence for the 
population of different quantum levels of the fluxon in the magnetic-field-induced potential 
well, which exists in the static state. Given that the model of this paper is classical, we 
must interpret our results differently. 



In the classical system under investigation here, we do not have intrinsic quantum energy 
levels of the macroscopic flux quantum. However, as in the experiments, distinct resonances 
between the fluxon oscillation and the ac field may appear when the ac frequency and 
the fluxon plasma frequency become harmonics of one another; specifically, when they are 

equal. The applied frequency u d = 0.2 used for the simulations of Figure|2l equals the fluxon 

f f) 

resonance frequency Up in the bias current vicinity of the secondary peak in the switching 
distribution for increasing ac amplitude. 

This is verified by Figure EJ in which we have displayed the relationship between the 
location of the secondary (the ac- induced) peak in bias current, and the predicted fluxon 
resonance frequency. To demonstrate that the multi-peak switching distribution in the 
classical limit is not limited to particular parameter ranges, we have chosen to show results 
for I = 20, T = 0.25, and fjdc — 2.5 ■ 10~ 7 . These parameters provide a critical fluxon bias 
current of rjc ~ 0.0699. Each distribution shown in Figure|3]is for a particular value of the 
driving frequency Ud and the corresponding ac amplitude has been optimized to produce 
multiple peaks of similar content in switching events. Figure El shows the predicted, Eq. (flj). 
relationship between the applied frequency and the location of the ac-induced peak as a 
solid curve. As can be seen, very reasonable agreement is found. We notice, however, 
that the peak location is consistently shifted slightly down in frequency (down in current) 
compared to the predicted linear resonance. This slight discrepancy is consistent for all 
system parameters we have tried (different lengths and magnetic fields). The reason for the 
variation is that the resonance leading to any switching event is not purely linear, since the 
fluxon must experience large oscillations in at least a cubic nonlinearity in order to escape the 
potential well. Thus, the resonance frequencies are depressed by third-order nonlinearities. 
This can be quantified from the following simple analysis. Let xq = Xq + A, where xq is a 
constant and A represents oscillatory motion of the fluxon around x . The non-relativistic, 
noiseless limit of Eq. then reads : 

A + —Vc^ s i n cos kA + — r]j/' cos kx sin kA = -^Vdc + ^Vac sin u d t — a A . (14) 

We assume a single-mode solution in the form 

A — a sm(uJdt + ijj) , (15) 
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which, when inserted into Eq. (JHJ), yields 



77 



sin kxr 



f] ac sm(u d t -ip) = a 



Vdc 



77 

2 



kr/P — !- — - cos kxo — 
ka 



sin Udt + aaujd cos Udt 



(16) 



(17) 



In the low dissipation limit, this results in the amplitude-dependent resonance frequency 

(18) 



(/) 



\ 



77, (/) Ji(feo) 



A 



Vc f) Uka) 



Notice that for small amplitude oscillations with \ka\ <C 1, Equations (fTTj) and (j!8)l coincide. 
For larger oscillation amplitudes a, however, we observe a decrease in both critical current 
and resonance frequency For multi-peaked switching distribution, as observed in Figure 01 
obtained for relatively low temperature 8, the dynamics must be such that the driven oscil- 
lations of the fluxon have an amplitude similar to the distance from the energetic minimum 
position to the corresponding energetic saddle position of the washboard potential; i.e., we 
can estimate the relevant amplitude of oscillations for multi-peaked distributions to be given 
by 



J (ka) cos ^ 



ka < 7i — 2 sin 

Vdc 



Vdc 



V { c f) J (ka) 



(/) ■ 



(19) 
(20) 



A polynomial expansion of the left-hand side of Eq. (|2T)|) for small ka yields the following 
approximate relationship between the fluxon oscillation amplitude and the bias current: 



ka 



\ 



4 


Vdc 


3 


JJ) 



(21) 



Inserting Eq. ()21|) into Eq. ()18j) provides an estimate of the anharmonic resonance frequency 
relevant for resonant multi-peaked switching distributions. This relationship, func- 
tion of f]dc, is shown in Figure El as a dashed curve. We clearly see the depression of the 
resonance and the very good agreement with the numerical simulations. We notice that 
fitting the critical bias r/P in the linear resonance curve shown in Eq. (fTTj) will lead to 
very good agreement between simulation results as well, but it is important to emphasize 
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that while Eqs. (jTHj) and (|21|) are approximate expressions they do not contain any fitting 
parameters. 

Further investigations of the ac-induced multi-peaked switching distributions and the 
relationship to the above perturbation analysis are shown in Figure 0] for a system with 
parameters / = 10, r\ dc = 10~ 6 , and r = |, |, |. We show here only the locations of 
the ac-induced switching distribution peaks as markers along with the corresponding linear 
resonance (Eq. (fTTj) . solid) and anharmonic resonance (Eqs. ()18|) and (}2"T]) . dashed). Again 
we find very good consistency between the simulated secondary resonant switching peaks 
and the approximate resonance curves. We see a slight discrepancy developing for increasing 
magnetic field T. This arises due to the strong sensitivity of the perturbation result to the 
predicted critical current as given by Eq. (jl(J|) . Figure shows the predicted critical bias as 
lines together with the corresponding values obtained from numerical simulations of Eq. 
for 9 = r] ac = 0. The agreement is good for most low to moderate values of T, but larger 
magnetic fields show that the predicted value overestimates the actual critical bias. However, 
even the very slight difference observed for I = 10 and F — | gives rise to the discrepancy 
observed in Figure 0] for T = |, and a fine tuning of 77c by less than 1% can make the 
agreement in Figure 0] be near perfect. 

The ac-induced resonant escape described above has been analyzed for situations when 
an applied ac frequency equals the fluxon oscillation resonance. However, one can equally 
well expect resonant escape when the applied frequency is a simple rational multiple of the 
resonance frequency, thereby exciting the fluxon oscillations. This is demonstrated by a few 
switching distributions shown in Figure |B1 for parameters I = 10 and T = 1. Here, q indi- 
cates the relationship between the driving frequency and the excited resonance q~ l uj d = u>rU, 
giving rise to the left-most (small bias) escape peak. As can be directly observed, both su- 
perharmonic and subharmonic resonances can result in ac-induced switching. Furthermore, 
we observe (for q = 2) that more than one resonance may be excited enough to provide more 
complex thermal switching distributions in some cases. 

Figure U\ summarizes an investigation of the effect of the sweep rate fjdc- Key parameters 
here are I — 10, Y — 1, u) d — 0.2125, and r\ ac = 10~ 3 . As should be expected, switching dis- 
tributions of different sweep rates differ such that lower sweep rates enhance early (low 77^) 
switching, simply because more switching attempts can be made for low bias values. Conse- 
quently, in the extreme slow sweep simulations (77^ = 10~ 8 ) we observe almost all switching 
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events near the (q = 1) resonance (rjdc ~ 0.3875). In the opposite limit (rjdc = 10 4 ), we 
observe nearly all switching events around (or beyond) the critical bias (rjdc ~ r)c ~ 0.41), 
and we also observe complex switching distributions for these very fast sweeps. In this case 
the sweep rate is too large for the fundamental resonance to develop within the time that 
rjdc provides fluxon resonances near the driving frequency, and the system consequently does 
not have significant opportunity to develop a large amplitude resonant state. Between the 
two extremes of Figure [?l we have distributions of sweep rates that provide distinctly sepa- 
rated switching peaks. We have, therefore, except for data shown in Figures [7| and |Hl chosen 
sweep rates that on the one hand provide multi-peaked distributions that are not artifacts of 
excessively fast sweep rates, and on the other hand provide simulation times that will yield 
statistically meaningful switching distributions with our available computer resources. We 
notice that the nearly single-peak distribution shown for rjdc = 10~ 8 in Figure 7 should not 
be interpreted as if multi-peaked distributions vanish for slow sweep rates. In fact, multi- 
peaked distributions for very slow current sweeps will become observable if the ac amplitude 
rj ac is slightly decreased. Likewise, the very fast sweep will reveal the resonant behavior 
if the ac amplitude rj ac is increased. Thus, we note that complex switching distributions 
may be easily obtained for many meaningful parameters in ac-driven Josephson systems. 
To emphasize this essential point, we have produced switching distributions for rjdc = 10~ 8 
for different values of r] ac . The results, shown in Figure |Hl clearly show how lowering the ac 
amplitude is compensating for the low sweep rate, thereby again illuminating the resonant 
switching behavior associated with the coupling to the applied ac field. 

In conclusion, we have thoroughly investigated the classical thermal sine-Gordon analog 
to the recently published experimental data from an ac-driven, magnetic-field-embedded 
annular Josephson junction with a single trapped fluxon. As a result of the applied ac per- 
turbations, we can observe resonant behavior of the trapped fluxon, which in turn gives rise 
to measurable escape signatures closely resembling experimental observations^ made in low 
temperature conditions. Simple considerations of the (anharmonic) resonances responsible 
for the anomalous switching distributions show very good and consistent agreement between 
simulation and analytic perturbation results. Since the external ac field creates resonances 
whose switching distribution signatures coincide with those of possible intrinsic quantum 
levels and quantum tunneling, we interpret the results as evidence that the identification of 
multi-peaked switching distributions cannot be unambiguously related to the identification 
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of quantum levels and quantum tunneling if the anomalous switching distribution is gener- 
ated with the application of an external ac field. This conclusion is further strengthened 
by recent experimental and theoretical worki^ demonstrating that multi-peaked switching 
distributions in small- area Josephson junctions can be induced by an ac perturbation for 
temperatures well above the quantum transition temperature similarly to what has been 
reported in the low temperature regime. 
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FIG. 1: (a) Sketch of the system under consideration: an annular Josephson junction, defined 
by the two superconductors Si and S2, with (normalized) bias rj = r]d c + ?7 ac sinu^i; (b) sketch of 
the trapped flux-quantum orientation and position relative to the external magnetic field H; (c) 
sketch of the magnetic-field-induced washboard potential felt by the flux quantum, which here is 
illustrated as a particle at position x = xq. 



13 




14 




FIG. 3: Simulated switching distributions p{r]dc) for I = 20, T = 0.25, rjdc = 2.5 • 10~ 7 , and 
9 = 10 -3 . Each distribution represents 5000 switching events. Top plot shows the linear resonance 
(solid curve), Eq. and the anharmonic resonance (dashed curve) of Eqs. Q18|) and (|21jl. Markers 
(•) show the corresponding frequency-bias relationships as obtained from the ac-induced resonant 
peaks of the displayed switching distributions. 
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FIG. 4: Simulated relationships between bias and ac-induced oscillation resonance for I = 10, 
-q dc = 10" 6 , and 9 = 2.5 • 1CT 3 . Each marker (T = 0.25 □; T = 0.5 O; r = 0-75 0) represents 
a distribution as shown in Figure with 5000 switching events. Solid curves represent the linear 
resonance, Eq. (|11|). and the dashed curves represent the anharmonic resonance of Eqs. ()18|) and 

(EH- 
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FIG. 5: Comparisons between the predicted critical bias (solid lines) from Eq. (|10|) and simulated 
values (I = 10 •; I = 20 Q) for 9 = Vac = 0. 
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P(Vdc) 




FIG. 6: Simulated ac-induced switching distributions for various harmonics (q) of the same fre- 
quency, uJd = Q ■ 0.225. Parameters are: I = 10, V = 1, fj^c = 10~ 6 , and = 2.5 • 10 -3 . Each 
distribution represents 5000 switching events. 
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FIG. 8: Simulated switching distributions p(% c ) for I = 10, F = 1, u>d = 0.2125, i]d c = 10~ 8 , 
and 6 = 2.5 • 10~ 3 for three different ac-amplitudes: r] ac = 10~ 3 (top), r/ ac = 7 • 10~ 3 (middle), 
rj ac = 5 • 10~ 4 (bottom). Each distribution represents at least 2000 switching events. 
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